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ON THE DISAPPEARANCE OF THE BROAD-LINE REGION 

IN LOW-LUMINOSITY AGNs 

Moshe Elitzur 1 and Luis C. Ho 2 
ABSTRACT 

The disk-wind scenario for the broad-line region (BLR) and toroidal obscu- 
ration in active galactic nuclei predicts the disappearance of the BLR at low 
luminosities. In accordance with the model predictions, data from a nearly com- 
plete sample of nearby AGNs show that the BLR disappears at luminosities lower 
than 5 x 10 39 (M/1O 7 M ) 2 / 3 ergs" 1 , where M is the black hole mass. The ra- 
diative efficiency of accretion onto the black hole is < 10~ 3 for these sources, 
indicating that their accretion is advection-dominated. 



Subject headings: accretion, accretion disks 
Seyfert — quasars: general 



galaxies: active — galaxies: 



INTRODUCTION 



According to the unification model for active galactic nuclei (AGNs) , the only difference 
between type 1 (broad-line) and type 2 (narrow-line) sources is the observer's orientation with 
respect to the toroidal dusty obscuration. Directions with clear sight of the central engine 
and the broad-line region (BLR) yield type 1 sources, while those blocked by the torus from 
direct view of the BLR result in typ e 2 objects, whe re the existence of the hidden BLR is 
revealed only in polarized light (e.g.. lAntonuccilll993l ). However, in spite of the consider able 
succes s of the unification scheme there is now clear evidence, recently summarized by iHo 
(120081 ). that the BLR is actually missing, and not just hidden, in many low- luminosity AGNs 
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(LLAGNs). These sources have been named "pure" flTranll200ll . 120031 ) or "true" flLaor!l2003r ) 
type 2 AGNs. Why does the BLR disappear in LLAGNs? 

In another important recent development, the obscuring torus is found to be a smooth 
contin uation of the BLR, not a separate entity. Primary evidence comes from the lSuganuma et al 
( 120061 ) discovery that t he BLR extends outw ard until the inner boundary of the dusty torus, 
in agreement with the iNetzer fc Laorl (19931) proposal t hat the BLR outer boundary is set 
by dust sublimation. Additionally, iRisaliti et al.l (120021 ) find that X-ray variability caused 
by passage of absorbing clouds across the line-of-sight has a continuous behavior across the 
time scales generated by the motions of dust-free BLR clouds close to the AGN and the 
more distant, dusty torus clouds. Together, these observations show that the BLR and the 
torus are simply the inner and outer regions, respectively, of a single, continuous cloud dis- 
tribution. Their different radiative signatures merely reflect the change in cloud composition 
across the dust sublimation radius Rd- The inner clouds are dust-free, and their gas is di- 
rectly exposed to the AGN ionizing continuum; therefore, the gas is atomic and ionized, 
producing the broad emission lines. The outer clouds are dusty, and so their gas is shielded 
from the ionizing radiation, and the atomic line emission is quenched. Instead, these clouds 
are molecular and dusty, obscuring the ultraviolet/optical emission from the inner regions 
and emitting infrared radiation. Within this framework, the BLR occupies r < R& while 
the torus is simply the region at r > Rd, and a more appropriate des ignation for it is the 
"toroidal obscuration region" (TOR; lElitzurl 120071 ; iNenkova et al.l 120081 ) . 



The BLR/TOR st ructure arises naturally in the disk-wind scenario, first proposed by 
Emmering et al.l (119921 ). In this model, the two classes of clouds simply correspond to dif- 
fere nt regions of a clumpy win d coming off the accretion disk rotating around the black hole 



isee lElitzur fc Shlosmanll2006l . and references therein). As the clouds rise away from the disk 
they expand and lose their column density, limiting the vertical scope of broad-line emission 
and dust obscuration and emission, resulting in a toroidal geometry for both the BLR and 
the TOR. Although a theory of clumpy disk winds in AGNs is far from full development, an 
immediate consequence of this scenario is the prediction that the TOR and BLR d i sappear 
at low bolometric luminosities (i.e., low accretion rates; Elitzur fc Shlosman 2006; Elitzur 



20081 ). The reason is that, as the mass accretion rate decreases, the mass outflow rate of a 
disk wind with fixed radial column decreases more slowly and thus cannot be sustained below 
a certain accretion limit. This unavoidable conclusion follows from simple considerations of 
mass conservation. However, in the absence of a complete theory for the outflow dynamics, 
the actual limit and its detailed dependence on the AGN parameters remain undetermined. 
Presently, the only practical approach is to attempt to extract clues on these unknowns from 
analysis of LLAGN data. This is the approach taken here. 
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2. The BLR Disappearance 



Denote by L45 = L/10 45 ergs -1 the AGN bolometric luminosity and by M7 = M/10 7 M Q 
the black hole mass. From mass conservation, disk outflow can support the column densities 
required by the BLR and TOR only when the luminosity obeys 



L 45 > CM. 



2/3 



X 



( lElitzur fc Shlosmanl 120061 ; lElitzurll2008l ). The derivation of this result is repeated below in 
§21 where we analyze in detail the assumptions and provide the expression for the numerical 
coefficient C in terms of system properties other than L and M. Replacing the mass with 
the Eddington luminosity LEdd = 1-3 x 10 45 M 7 ergs _1 , produces the luminosity lower bound 
in terms of Eddington ratio 

-^Edd 



L45 > C<>, 



(2) 



where C2 = 0.59 C 3 . In the absence of a detailed theoretical model for the disk wind, our 
ignorance about the dynamics is contained in the numerical coefficient C. The only way to 
determine whether C may contain additional, indirect dependence on L or M is to compare 
the bounds in eqs. [1] and [2] with observations. 

We use dat a selected from the Palomar spectroscopic survey of 486 nearby galaxies 
f lHo et dJE997aB l. the most sensitive and complete sample of LLAGNs available. A critical 



(Ho 


2009alJbl: 


Ho et al. 


2009) 



all galaxies classified as Seyferts, low-ioni zati o n nuc l ear emi ssion-line regions (LINERs), and 
transition objects, which, as argued in 155 J2008I . l2009ah . are variants of LLAGNs that 



define a sequence of varying accretion rates. We also include absorption-line nuclei, objects 
with undetected optical line emission that are good candidates of hosting inactive black 
holes. High-resolution, hard-X-ray measurements provide the most robust measure of the 
accretion luminosity in AGNs, especially for weak sources. Of the 277 Palomar galaxies 
meeting our classificatio n criteria, 166 (60%) have s uitable nuc lear X-ray (2-10 keV) fluxes 



or upper limits thereof (IHol l2009al ) . As discussed in iHol (j2009al ). these X-ray measurements 
can be converted to bolometric luminosities with an accuracy of ~ 0.3 dex. Importantly, all 
but one of the galaxies has reliable measurements of the central stellar velocity dispersion 
( iHo et al.ll2009l ). from which the black hole mass can be estimated to within ~ 0.3 dex using 



the well-established re lation between black hole mass and bulge stellar velocity dispersion 



(ITremaine et al.l 120021 ). 



Figure 1 shows the distributions of data points in the L — M and L — L/Lsdd planes. 
The Palomar LLAGNs are coded according to spectral class and presence or absence of 
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broad Ha emission ( iHo et al.lll997bl ). For comparison, we also include the sample of low- 
redshift luminous Seyfert 1 nuclei and quasars studied by lGreene &: Hoi (120071 ). It is apparent 
that, depending on either the black hole mass or Eddington ratio, type 1 sources cease to 
exist below a certain luminosity (The only exception is the LINER 1 galaxy NGC 2787, 
whose luminosity of 10 39 ' 98 ergs -1 is 0.45 dex below the drawn boundary line.) We interpret 
this critical luminosity to be the threshold below which the BLR disappears in the disk- 
wind model. The solid line in each panel represents the bound with the slope taken from 
either eq. [T] or [2J as appropriate. The line intercepts were adjusted to coincide with the 
observed type 1— type 2 boundary, yielding C = 4.7 x 10 -6 , i.e., the BLR existence requires 
erg s" 



L> 5 x 10 39 — 



3. The Disk- Wind Constraint 

Consider an outflow coming off an annular segment of a disk. Denote by n(R) the density 
and v z (R) the vertical component of outflow velocity at axial radius R on the disk surface. 
Then the outflow mass loss rate is M out = lism J nv z RdR, where m is the mean atomic mass. 
Since the dust sublimation radius Rd sets a characteristic scale for both the BLR and TOR, 
it is convenient to introduce the scaled radius y = R/Rd and the dimensionless velocity 
profile u = v z (R)/v z (R^). Introduce also N R = J n(R)dR, the overall column density along 
a radius vector on the disk surface. We can describe the radial density variation by the 
dimensionless profile r\ = n(R)Rd/N R , normalized from J r](y)dy = 1. Then 

M out = 27rmN R R d v z (R d )I , (3) 

where 

I = J v{y)u{y)ydy 

is a numerical factor of order unity. Written in terms of the radial column density Nr, this 
basic expression provides the mass outflow rate of any disk wind, whether smooth or clumpy; 
in the latter case v z is the cloud liftoff velocity and the profile rj describes the variation of 
number of clouds per unit radial length. 



1 /2 

The dust sublimation radius in AGNs is Rd — 0.4L 4 5 pc (INenkova et al.l 120081 ). The 
initial outflow velocity, v z , is roughly the dispersion velocity of material in the disk — an 
outflow is established when the ordered motion velocity becomes comparable to that of the 
local random motions. We now make the assumption that this velocity is some fraction fi 
of the local Keplerian velocity, v z (Rd) = fi(GM/ Rd) 1 ^ 2 . Maser observations of the nuclear 
disk in NGC 3079 show that the velocity dispersion is ~ 1 4 kms" 1 in a small re gion of 



strong emission where the Keplerian velocity is 110 kms 1 (IKondratko et al.ll2005l ). This 
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observation suggests that fi is probably of order ~ 0.1. With this assumption, the mass 
outflow rate becomes 

Mont = 0.07 fJN 22 M^ /2 L\i 4 M G yi~\ (4) 

where N 22 = Nr/10 22 cm -2 . This expression applies both for the BLR and TOR portions of 
the disk wind. 

In steady state, M out cannot exceed M acc , the mass accretion rate into the BLR/TOR 
region. A fraction f 2 of this accreted mass finds its way to the black hole and is converted to 
the AGN bolometric luminosity with radiative efficiency rj so that L = f 2 f]M acc c 2 . Therefore, 
M acc = 0.02^45/(^/2) Moyr" 1 and 



M, 



out 



Mr 



1/2 



L 



3/1 

15 



(5) 



where e = 3.5rjfif 2 IN 22 . Steady-state mass conservation requires M out /M acc < 1, yielding 
eq. Q] with C = e 4/3 . The observational result C = 4.7 x 10~ 6 implies T]fif 2 IN 22 = 3 x 10" 5 . 
From photoionization mo deling of broa d-line emission, the BLR requires a minimal column 
density of ~ 10 22 cm" 2 (INetzeri Il990h . and so iV 2 2 > 1. We can reasonably take I ~ 1 
and /1 ~ 0.1, yielding f 2 r] ~ 3 x 10 -4 ; since the factors f 2 and i] do not enter separately, 
only their product is constrained. A rough estimate of f 2 can be obtained from self-similar 
hydromagnetic disk-wind models. The wind applies a back torque on the underlying disk 
because each magnetic field line can be considered rigid from its footpoint R up to the 
Alfven radius R\. Angular momentum conservation implies that 7 = (R/R^) 2 is roughly 
the ratio of the mass fluxes for outflow and accretion, or 7 ~ M out (i?)/M acc (i?). This ratio 



i s found to be either inde pendent of R (jEmmering et al.lll992l ) or weakly dependent on it 
( jPelletier fc Pudritzlll992l ). with a value of 7 ~ 0.1 — 0.3. The fraction of M acc that reaches 
the black hole is 1 — 7 = f 2 , therefore f 2 > 0.7 and rj ~ 4 x 10~ 4 . Such a low radiative 



Naravan 


2002; 


Yuan 


2007) 



"standard" accretion efficiency 77 ~ 0.1 would require virtually the entire mass accreted into 
the BLR/TOR to be carried away by the wind, with only a fraction f 2 ~ 3 x 10~ 3 reaching 
the black hole. This seems unlikely. 



4. Discussion 

Quenching of the BLR/TOR at low accretion rates is an unavoidable consequence of 
the disk-wind scenario. The reason is that the mass outflow rate is set by the radial column 
density N R and by the scales of the local radius and Keplerian velocity (eq. |3]). Since the 
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BLR/TOR boundary is set by a fixed value of the radiative flux, the radial scale increases 
with luminosity as L l l 2 and the velocity scale decreases as L~ 1//4 . Keeping Nr constant, as 
appropriate for both the BLR and TOR, the mass outflow rate varies with L as L 1 / 4 , more 
slowly than the linear variation of the accretion rate. Disk-wind mass conservation implies 
that the minimal column required by an observable BLR cannot be sustained below a certain 
luminosity. The data verify that this is indeed the case and indicate that the transition to 
"true" type 2 AGNs occurs in the regime of radiatively inefficient accretion. This adds 
sup port to independent indicators that LLAGNs accrete in the advection-dominated mode 
(|Holl2009al . and references therein). 



Bounds similar to eq.[T]were derived by lNicastrol (120001 ) and lLaorl (120031 ). Both concluded 
independently that the BLR should disappear when the luminosity drops below some L min = 
CM a , where C and a are constants. Replacing M with the Eddington ratio, this condition 
produces significantly different results depending on the value of a: when a < 1 it implies 
a lower limit L > C 2 (L^dd/ L) 1 ^ , where C% is derived from C, but when a > 1 the result 
is an upper bound L < G% (L / LEddl"" 1 j i n between, a = 1 implies that the BLR exists only 
above some fixed Eddington ratioa Nicas tro suggest s a dis k outflow origin for the BLR, as 
here, and uses the disk/corona model of IWitt et al.l (119971 ) to derive a = 7/8, practically 
indistinguishable from the a = 2/3 found here. However, Nicastro's bound translates into 
L > C2(L-Edd/L) 7 , a much steeper boundary than indicated by the data in the L — L/L-^dd 
plane (Figure 1). Laor suggests the existence of an upper limit to BLR velocities, which 
leads to a = 2. Because a > 1 in this case, it yields L < C 2 (L/L Edd ) 2 . This bound is 
orthogonal to the one plotted in panel (b) of Figure 1 and all AGN with BLR should lie 
below it, a prediction that seems in clear conflict with the data. As these examples show, 
examining the data in the L — L/L-^dd plane amplifies the differences between models; while 
eq. [2] describes adequately the data in Figure 1(b), the other models do not. 

The derivation of the BLR/TOR disapp earance presented here is more general than 
the original one in lElitzur fc Shlosmanl (120061 ). Equation H] holds irrespective of the wind 
dumpiness and applies equally to the BLR and TOR portions of the disk outflow, with 
N 2 2 the radial column density through each pertinent segment. In both regions, N 2 2 ~ 1 
can be considered a reasonable lower limit for the BLR to generate line emission and the 
TOR to produce substantial obscuration. The only significant difference between the two 
cases enters from the factor I (eq. [3]): denoting the inner and outer radii of the disk wind 
by -R; n and -R out , respectively, the integration in /blr extends from i? in /i?d to 1, while in 
/tor it is from 1 to R out /Rd- We can reasonably assume Keplerian behavior for the velocity 



1 As a increases, the bound is rotating clockwise in the logL — log(L/LEdd) plane, passing through vertical 
at a = 1, with BLR sources located ahead of the rotating boundary. 
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profile, i.e., u 



y 



-1/2 



and parametrize the density profile with a power law, rj oc y~ 



so that the integrand is proportional to y 1 ^ 2 q . Depending on q, different relations exist 
between /blr and /tor- The range f < q < 3/2 yields Itor/^blr ~ 1, Q < 1 produces 
/tor/W « (^out/^d) 172 > 1, while g > 3/2 results in / TO r/W « (iW^d) 1/2 < 1. 
The larger is /, the earlier does the quenching of the outflow occur as the luminosity is 
decreasing. Therefore, depending on the steepness of the radial density profile, the BLR 
and TOR may disappear either together, when Jblr ~ -^tor, or i n sequence: in AGNs 
with steep profiles (q > 3/2; /rlr > -^tor) the BLR will disappear first; in shallow ones 
(q < 1; /tor > /rt,r) the TOR. Recent clu mpy torus modeling of extended IR data sets 
( IMor et al.ll2009t iRamos Almeida et al.ll2009l ) suggests that q can range anywhere from to 
3 in different sources, indicating that no single rule applies to all AGNs — sources with and 
without TOR can exist on either side of the BLR boundary. In any given AGN, the TOR 
disappearance can be established from the lack of IR emission at a level commensurate with 
the bolometric luminos i ty. Th e absence of a TOR in M87 was demonstrated conclusively by 



Whysong fc Antonuccil (120041). who place d stringent limits on the thermal IR emission, and 



further solidified by lPerlman et al.l (120071 ). who found only a trace of thermal emission that 
is much weaker than expected from an AGN torus and that can be attributed to neighboring 
dust. Detailed studies of the IR emission from sources close to the BLR boundary established 
in Figure 1 can be a powerful tool for probing the density distribution around the black hole 
in LLAGNs and clarifying the nature of the BLR/TOR disk outflow. 
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Fig. 1. — Distribution of (a) black masses, M, and (b) Eddington ratios, L/L Edd , vs. bolo- 
metric luminosity, L, for objects separated by spectral classification. The bolometric luminos- 
ity is based on the 2-10 keV X-ray luminosity, as described in Ho (2009a). The symbols are 
identified in the legend. Solid and open symbols indicate type 1 and type 2 sources, respec- 
tively. "S" = Seyferts, "L" = LINERs, "T" = transition objects, and "A" = absorption-line 
nuclei. Short line segments denote upper limits. The objects marked as "QSO" refer to the 
sample of high-luminosity Seyfert 1 nuclei and quasars studied by Greene & Ho (2007). The 
solid lines are logL = 35 + |logM in panel (a) and logL = 28.8 — 2 log(L/LEdd) in (b), 
corresponding to the theoretical predictions given in eq. (1) and eq. (2), respectively. 



